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Abstract 

q | We consider the classical theory of a gravitational field with spin 2 and 

(Pauli-Fierz) mass m in flat spacetime, coupled to electromagnetism and 

|h . point particles. We establish the law of light propagation and calculate the 

amount of deflection in the background of a spherically symmetric grav- 
itational field. As m tends to zero, the deflection is shown to converge 
k> \ to 3/4 of the value predicted by the massless theory (linearized General 

;_j ' Relativity), even though the spherically symmetric solution of the gravi- 

tational field equations has no regular limit. This confirms an old argu- 
ment of van Dam and Veltman on a purely classical level, but also shows 
its subtle nature. 

Notation and Conventions 

We work in flat Minkowski space with metric r\ = diag(l, —1,-1,-1). Space- 
time indices, running from to 3 are greek, whereas space indices, running 
from 1 to 3 are latin. The Minkowski inner product is denoted by rj{p, q) — 
P ■ q = p^q^, andp ■ p = p 2 . Partial differentiation with respect to x^ is either 
written as d^ or sometimes simply as , \x. Throughout we use units in which 
the velocity of light is set to 1 . 



Introduction 

According to an old and well known argument of van Dam and Veltman [|2| |3| , 
which we will sketch below, the theory of massive spin 2 fields in flat 
Minkowski spacetime does not approximate the strictly massless theory (lin- 
earized General Relativity (GR)), in the limit as the mass tends to zero. (See 
also [§] and the more comprehensive account in (j} .) This means that there 
exist corresponding observables in both theories which distinguish the mass- 
less limit of the massive theory from linearized GR. One such observable, so 
it is claimed, is the amount of deflection of a light ray in the background of 
a, say, rotationally symmetric, static gravitational field. More precisely, in the 
limit the mass tends to zero, the deflection angle predicted by the massive 
theory tends to 3 /4 of the angle predicted by linearized GR. 

As we shall see below, the argument given by van Dam and Veltman is 
entirely based on the properties of free propagators, whose structure is deter- 
mined by Poincare invariance. The aim in this paper is to improve our under- 
standing of the classical aspects of this limit, which is known to be non trivial 
for several reasons, but has not been explored in all details in the mentioned 
references. In particular, we wish to understand how certain observables can 
have a smooth limit as the mass tends to zero, even though they refer to solu- 
tions which diverge in that limit. The deflection of light is an example of such 
a case, whose derivation from first principles we shall consider in this paper. 

Let us now turn to the argument proper of van Dam and Veltman: In mo- 
mentum space, the free propagators are given by 

P™ aP (p) = 5(WW» + *M/»^J-§W«/» massiveca se (1) 

^ p z — m z + ie 

P^(P) = Kw^ + ^)-K^ mass .iesscase (2) 

p A + ie 

where n^ is the induced Riemannian metric on the mass-hyperboloid p 2 = 
m 2 in momentum space: 

■Kpvip) ■= - JL -^- - V/MU ■ (3) 

We now consider two systems with (Fourier transformed) energy-momentum 
tensors T^ v and T'. If we assume these tensors to be conserved, = 
p^Tftvip) = p tl T'^(p), the one-graviton interaction takes the form (we write 
T = T£ etc.) 

.*~W*5-M !•*(-»> - « ^M^rtzinggHj (4) 

p M p z — m + ie 



in the massive case, and 

k T m {p)P t _ iya p{p)T (-p) = k c 3-— (5) 

p z + ie 

in the massless case, with gravitational constants k and k respectively. 
For slowly-moving massive objects the leading component of the energy- 
momentum tensor is Too and we get |kT 00 (p)Tq (-p) in the massive and 
^k T 00 (p)Tq (— p) in the massless case. To identify k and k one requires that 

both cases lead to Newton's law of attraction with the same Newtonian con- 
stant G (taking the limit m — ► in the massive case). This gives 



3,. 



«0 



VVI^G. (6) 



On the other hand, if we consider the interaction of light and slowly mov- 
ing matter, there are no trace-terms in (§) and (g) due to the tracelessness 
of the electromagnetic energy momentum tensor. But now @ implies that 
the massive theory, in the limit m — ► 0, leads to an interaction of light and 
matter which is weaker by a factor 3/4 as compared to the massless theory. 
Accordingly, in that limit, light deflection comes out smaller by the factor 3/4 
as compared to the massless theory, i.e., linearized General Relativity. This 
is often taken as sufficient justification to rigorously argue that present-day 
observations exclude a finite graviton mass. 

Here one should stress that the flatness of spacetime enters in an essential 
way. In fact, it was recently argued that there is no discontinuous behaviour 
at to = if the cosmological constant A is non-zero. For negative A (anti 
de Sitter spacetime) a smooth limit of propagator residues was shown in [||] . 
For positive A (de Sitter spacetime) the situation is again different, since there 
is no quantum theory of spin 2 fields in the mass range < m 2 < | A obeying 
unitarity and certain locality requirements [§|, 0] . Here a discontinuity occurs 
atm 2 = | A. 
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1 Massive spin- 2 fields in Minkowski space 

We briefly review the theory of massive spin-2 fields in flat spacetime. It is 
represented by a symmetric tensor field h^. The space of such fields still 
represents the Poincare group in a highly reducible fashion. A projection into 



an irreducible subspace of pure spin-2 and mass m is given by 

(□ + m 2 )V = (7) 

0"V = (8) 

h* = (9) 

1 . 1 Static rotationally symmetric solution 

We seek a solution of (0,§J9l) which is static, i.e. d a h^ = and rotationally 
symmetric in the following sense: 

D{R)«D(R)l h afj {D{R) x a x°) = h^(x x ) WR e 50(3) , (10) 

where, in a 1+3 - split matrix notation, 

For physical reasons (finite energy around spatial infinity) we also require the 
asymptotic fall-off condition lim r _,oo h^ v (x) = 0, where r is the 3-dimensional 
radius. The 00-component of ^ has the one-parameter (b) family of solutions 

hoo(x) = -bf(r) := -6 eXp( ~ mr) , (12) 

r 

where the minus sign is introduced for later convenience. The 0-component 
of (||) and ( |To| ) imply h 0l = ax l /r 3 , which contradicts (Q) unless a = 0; hence 
hoi = 0. To determine the spatial components, we first remark that any rota- 
tionally symmetric two-tensor in space is of the form 

hij(Z) = fi(r)Sij + f 2 (r)°^ . (13) 

Equation ^ now reduces to two coupled ODEs for f\ and / 2 , which may be 
decoupled by introducing the new function /i := 3/i + f 2 . One obtains 

(A-m 2 )/i = 0, (14) 

(A-m a -£)/ 2 = 0, (15) 

which, under the given fall-off conditions, have the unique 2-parameter set 
of solutions 

fi{r) = Ci , (16) 

r 



exp(— mr) f 3 

mr ' (mr)' 



/a(r ) = c 2 ^^( 1 + _^ + _^_). (17 ) 



Condition Q and @ now imply c 2 = —\ci and c\ = -b respectively, thereby 
projecting out a unique one-parameter set of solutions, which, in terms of the 



function / defined in ( 12), can be written in the simple form 



h^x) = -b 



f(r) T 

|(^/(r) - ^a»fi,/(r)) 



(18) 



We note that the didjf - part is of a form that would be pure gauge in the 
massless theory. However, there is no gauge freedom in the massive theory 
and the prefactor, mT 2 , causes this term to diverge as m — > 0. 



A slightly more geometric way to write the solution ( |18[ ) is as follows: Let 
n — x/r be the radial unit vector; we define the spatial projection tensors 
pij := riirij in radial direction and m := % — p %J in the orthogonal direction, 
tangential to the spheres of constant r. In terms of these, the spatial part of 



( |18| ) takes the form 

hij(S) 



bf(r) 



(mr)' 



- (jij Lpij ) 



which clearly separates the trace part oc r, which stays finite for m 
the trace free part a r - 2/5, which diverges as m tends to zero. 



(19) 



0, and 



1 . 2 Lagrangian formulation and matter couplings 



The previous equations (|7|J8||9j) are equivalent to the Euler-Lagrange- 
equations of the following action: 



^S — 4 



d 4 x 



K v , x h^ x - 2h* xX x h llu > v + 2W v iU h tl 



m 2 {K v h^ - h 2 ) 



h^h'i* 



(20) 



The mass term is sometimes called the Pauli-Fierz term. A straightforward 
calculation of the variational derivative E^ := SSg/Sh^ (keeping in mind 
the symmetry h^ — h Ufl ) shows that the conditions E**" = 0, d^E^ = 0, and 
E£ = indeed imply (0,§,||) and vice versa, given that m ^ 0. 

The coupling to some specific form of matter is described by an interac- 
tion term 



S int = -^j 'd*xh llv T>™, 



(21) 



where T^ v is the energy-momentum tensor of the matter. The possibility of 
an /iT-coupling is excluded by an argument given in section (1.4). Taking the 



divergence of 6(S g + S [ s nt )/5h f ^ l/ we get 

m 2 d ll (h^ - rf v h) = -Kd^T^ . (22) 

Let us assume energy- momentum conservation for the matter field 

d^ u = . (23) 

Then d^W v = h' u , which, upon insertion into the Euler-Lagrange equations 
and taking the trace, implies that the traces of h^ v and T^ v must be propor- 
tional 

h=^T. (24) 
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Using (E2fi and (24) in the Euler-Lagrange equations they read 



-(□ + m 2 ) V = K (V - l VlJlv T) - -^dpdvT . (25) 

This leads indeed to the propagator ([j} and establishes the background for 
the arguments given in the introduction. 

1 .3 Coupling to a point particle 

The action of a free point particle of mass m in Minkowski space is given by 
— m times its arc-length. Choosing any parameter A to parametrize its world- 
line, z M (A), we have 



S p = -mo / d\ Wr^ 



dz^(X) dz u (X) 



dX dX 
Its energy-momentum tensor is 



~„„, n f , r ./A\ , , nn dz^ (t) dz" '(t) 

T£ v {x) =moldT 5^ {x - z{t)) ^ ; ^' 



(26) 



(27) 



where the parameter r is the arc-length with respect to the Minkowski metric 
r). Since this expression is not reparametrisation invariant, we are not free to 
specify it otherwise. However, we can play the following trick, which turns 
out to be of central importance: Consider a new metric on Minkowski space, 
given by 

fiV = V^ + Kh nv (28) 

The arc length w.r.t. g is called s. We consistently neglect terms of higher 



than linear order in k and may therefore replace r by s in (27), since this term 



already gets multiplied by k in (pip. In ( |26| ) we choose A = s, replace r\ by 
g - nli, and expand to linear order on k. The K-term then just cancels the 
interaction term ( pl| ) and we get 

S p + ff* = -mo / ds ^ g ^( z (s)) ^2 S) ^^ = ~mo J ds , (29) 

which says that the gravitationally interacting particle moves on geodesies in 
the metric <?. 



In order to determine the still unknown constants k and b in g^ v (see (181) 
we compare now this equation of motion with the Newton equation. To first 
order in k this reads, writing a dot for differentiation with respect to s: 

& = %rT (d v h a p - d a h vfi - d p K a ) z a z? . (30) 

In a Newtonian approximation, where we neglect terms of second and higher 
powers in the spatial velocities dz k /ds, this reduces to 

z° = - K h 00 , k z°z k , (31) 

z k = -ffcoo,*(*°) 2 - ®2) 

Introducing the Minkowski time coordinate t := z° we can rewrite the sec- 
ond equation, using the first one and again neglecting quadratic terms in the 
spatial velocities, as : 

d 2 z 



dt 2 



$Vh 00 . (33) 
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Thus, recalling the fall- off condition for /i M „, we get the identification 

^oo = *, (34) 



where $ is the Newtonian potential. Applied to our solution ( |12| ) we can now 
determine the product of the constants k and b to be 

nb = 2GM, (35) 

where G is the Newtonian gravitational constant and M is the central mass, 
i.e., the source of the gravitational field. With this identification our metric 
coefficients g^ v are now entirely determined. 

At this point we may also determine the coupling k separately. For this we 
return to the field equation (EJ3), whose 00-component for a static source and 
in the limit m — * reduces to 

A/too = §« T 00 . (36) 



In view of (|34|), a comparison with the Newtonian equation A$ = AirGp (here 
p = Too) leads to 

k = Vl27rG. (37) 

Finally, we also remark on the fact that instead of (||) one often finds the 
following alternative expression for the action of a free particle (e.g. in J9|): 

, f „ dz^(X) dz u (\) 

S p = -±m / dXrj^— ^-1—^-1. (38) 

Like ( ^ ) it also leads to straight lines in Minkowski space, but in addition 
also to the condition that A is (a constant multiple of) r, the arc length mea- 
sured with r\. Since one is only interested in the spacetime paths and not in 



their parametrisation, one may use either (26) or (B8J) to determine the free 



paths. However, we wish to point out that in the presence of an interaction 



of the form (21,270 the latter choice becomes inconsistent. The reason be- 



ing the following: Adding the "free action" (38) with A = r to the interaction 



(21) with the energy momentum tensor being given by ( 27 ) , we obtain an ex- 
pression like (^), where r) is replaced by g and the parameter is still r. But 
as a result of the fixed parametrisation the Euler- Lagrange equations now im- 
ply h^v dz ^ dz d i = 0> which in general will have no solutions. Hence we 
maintain that the action for the particle should be written in the square-root 
form @. 

1 .4 Coupling to the electromagnetic field 

The action for the free electromagnetic field in Minkowski space is given by 

S cm = -^Jd i x rf^ (iv F aP F^ , (39) 

where F^ u = d^A v - d v A^. The corresponding energy-momentum tensor is 
TZ = -F^F\ + \rT F a0 F al3 , (40) 



which we couple to the gravitational field according to ( |2l[ ) . A priori, there 
is no reason why different types of matter must couple with the same con- 
stant k. However, if we choose different constants for different matter types, 
say ki for type 1 and k 2 for type 2, then the ratio of gravitational to inertial 
mass of type 1 would differ from the corresponding ratio of type 2 by a factor 
K1/K2 and hence violate the weak equivalence principle, according to which 
the ratio of gravitational to inertial mass is the same universal constant for all 
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types of matter. Hence, in order to conform with this principle, we choose the 
same k for point particles and electromagnetism. This is also the reason why 



we excluded a hT- coupling in (21 1. 



Now we can play the same trick as for the point particle and write |g] 

S cm + S^ = -\Jd*xrf*rf v F afi F l „,-zJ#xh llv T' a ' 

= -\jd A x V- det{ S } g a V F^F^ + 0(« 2 ) , (41) 

where g^" and det{g} are the inverse and determinant of the matrix g^ LV . To 
first order in k we have g 1 *" = rf v — kW v and yj — det{.g} = 1 + |/i, where 
h = M and indices on h are always raised and lowered with g. 



We are now in the same situation as for the point particle, since Q4l[ ) is the 
action for the electromagnetic field on a spacetime with background metric 
g. Without sources its equations of motion are 

V^A] = ° ^ (V^ifs^ Fj) = , (42) 

which, in an eikonal approximation, imply that light rays are lightlike 
geodesies in the metric g. Note that our decision to couple all types of matter 
with the universal coupling constants has the effect that light rays and parti- 
cle paths are geodesies in the same metric. 



2 Deflection of light 

In this section we calculate the deflection of a light ray by the spherically sym- 
metric gravitational field (^8|). We are interested in the amount of deflection 
as m tends to zero. A priori this limit is a precarious one since the metric coef- 
ficients diverge as m — > (see (p|)). However, we now show that the deflection 
has a finite limit. To this end we split ( |T^ ) in a finite part, 

^w-'CT i4i)' (43) 

which has a continuous limit as m — > 0, and a diverging part, 

h%(x) = bm- 2 d i d j f{r). (44) 

Since we strictly stay within the linear theory and, accordingly, keep only 
terms linear in k, the contributions of h { and h d to the deflection add lin- 
early, so that we can consider them separately. Expanding /(r) = cxp(-mr)/r 



in powers of m, we explicitly checked that no power smaller than 2 of m 
in Q) contributes to the deflection, so that in the limit as m — ► the di- 
verging part does not contribute at all. In fact, one can argue that for each 
m the whole diverging part does not contribute to light deflection. The ar- 
gument is as follows: Recall that under a spatial coordinate transformation 
x l h^ x' 1 — x l + nk l (x j ) the spatial components of the metric g^ change - 
to first order in k - according to g i3 i-» g tj — K(kij + kj t i). Hence, choosing 
fc,(x J ) = hbm~ 2 dif(r), we can remove h d from g = t] + K{h { + h d ). The re- 
sulting coefficients g MV = i]^ + nh 1 ^ are then understood with respect to the 

coordinate system x n (we drop the dash hereafter). Since fc' falls of as r — ► oo, 
the new coordinate axes asymptotically approach the old ones. In particular, 
they are again asymptotically orthogonal. Hence we may calculate the light 
deflection 8{m) as usual, using the metric g i = i] + kH { with the identification 
(35). S(m) is a continuous function of m through the continuous dependence 



of the light deflection on the background metric. Hence, in order to obtain 
6(0), we may calculate the light deflection for the limit metric 

2GM qt 

r -5,,, (I - 



lim n sjL,0?) = 

m— »0 ^ 



l 



Ghr 



(45) 



But instead of performing an explicit calculation, we merely need to compare 



Q45P with the Schwarzschild metric in General Relativity. In isotropic coordi- 
nates the latter reads 



Schw 



(X) = 



1 + 



1- 



TTJT 



■2<;m 



T 







[1 

fl- 



GM ] 4 
2r J 



2GM \ 



C(G 2 ). 



(46) 



(47) 



As is well known, in leading (linear) order of G, light deflection receives equal 
contributions from the time-time and space-space parts: 

2GM 2GM AGM 



srEinst 



(48) 



q q q 

where q is the impact parameter. Applied to our metric this finally leads to 
the van Dam-Veltman value for the light deflection in the zero-mass limit: 

2GM GM 3GM 3 



lim S(m) 



cEinst 



q 



(49) 
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